Role of three-nucleon forces and many-body processes in nuclear pairing 
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We present microscopic valence-shell calculations of pairing gaps in the calcium isotopes, focus- 
ing on the role of three-nucleon (3N) forces and many-body processes. In most cases, we find a 
reduction in pairing strength when the leading chiral 3N forces are included, compared to results 
with low-momentum two-nucleon (NN) interactions only. This is in agreement with a recent energy 
density functional study. At the NN level, calculations that include particle-particle and hole-hole 
ladder contributions lead to smaller pairing gaps compared with experiment. When particle-hole 
contributions as well as the normal-ordered one- and two-body parts of 3N forces are consistently 
included to third order, we find reasonable agreement with experimental three-point mass differ- 
ences. This highlights the important role of 3N forces and many-body processes for pairing in nuclei. 
Finally, we relate pairing gaps to the evolution of nuclear structure in neutron-rich calcium isotopes 
and study the predictions for the 2"*" excitation energies, in particular for ^''Ca. 

PACS numbers: 21.10.-k, 21.30.-x, 21.60.Cs, 27.40.+Z 
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I. INTRODUCTION 

Pairingis a basic component in theoretical descriptions 
of nuclei [lll21- Pairing correlations, experimentally man- 
ifested in odd-even mass staggering, can be studied from 
three-point mass differences, which offer the best com- 
promise between capturing theoretical pairing gaps and 
excluding contributions from neighboring nuclei. The 
three-point mass difference for neutrons A„ is given by 



A(3) = 



-1) 



N 



-[B{N + 1,Z) + B{N -1,Z)- 2B{N, Z)] 



where B{N, Z) denotes the (negative) ground-state en- 
ergy of a nucleus with N neutrons and Z protons. 
This allows comparison between different theoretical ap- 
proaches and with the experimental A„ values obtained 
directly from mass measurements. Indeed, experimental 
investigations of exotic nuclei are advancing our under- 
standing of pairing towards the limits of existence |3|. 
Moreover, since a strong peak in the three-point mass 
difference indicates a particularly well-bound or magic 
configuration, it provides information for understanding 
the evolution of shell structure. 

Pairing gaps and global properties of medium-mass to 
heavy nuclei are best accessed using nuclear energy den- 
sity functional (EDF) methods 0|. The present EDFs 
are very successful, but are mainly constrained by fits 
to nuclei close to stability and therefore can lead to un- 
controlled uncertainties in experimentally unexplored re- 
gions, especially towards neutron-rich extremes [5|. As a 
result, recent efforts have focused on the construction and 
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constraint of EDFs based on microscopic two- (NN) and 
three-nucleon (3N) forces |6l-[ll|. Important steps have 
also been made in a framework focused on a microscopic 
treatment of pairing developed from NN and 3N interac- 
tions, while using existing empirical Skyrme functionals 
for the remaining part of the EDF |I2h15| . 

Pairing correlations and their influence on the proper- 
ties of medium-mass nuclei have also been studied within 
the nuclear shell model [l6l - t20| . To study nuclei towards 
the limits of existence, we have recently developed mi- 
croscopic calculations of valence-shell Hamiltonians using 
many-body perturbation theory (MBPT) based on chiral 
NN and 3N forces fit only to few-nucleon systems |2l1425| . 
In this approach, it has been shown that 3N forces play a 
key role for the evolution to neutron-rich and proton-rich 
nuclei [26| : for the location of the neutron dripline in the 
oxygen isotopes [2l|, [23| , for the ground-state properties 
and spectra along the proton-rich Z — ^ and Z — 2Q 
isotones [23], and for the A^ = 28 magic number in cal- 
cium [11] and the masses of ^^'^^Ca [2^]. The impact of 
3N forces on neutron-rich nuclei has also been demon- 
strated in coupled-cluster calculations with phenomeno- 
logical 3N forces |27l . l28ll. in the in-medium similarity 
renormalizati on g roup |29l . ISOj , and using Green's func- 



tion methods [31] 
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An important question regarding pairing in nuclei 
is the contribution of many- body processes, such as 
medium polarization i32]435|. When Coulomb effects 
were explicitly included in EDF calculations jl3,], there 
was a surprising agreement with experimental data al- 
though higher NN partial waves fSS!], 3N forces, and the 
contributions from many-body processes were neglected. 
When the contribution of 3N forces to pairing interac- 
tions was recently explored for the first time in EDF 
calculations [15|, this picture, however, changed signifi- 
cantly. It was found that 3N forces have a pronounced ef- 
fect, reducing the pairing gaps systematically by ^ 30%. 



As a result, the consistent inclusion of many-body pro- 
cesses and whether they can account for the missing pair- 
ing strength remain open problems in EDF calculations. 
In this work, we address this question in microscopic 
valence-shell calculations of neutron-rich calcium iso- 
topes that include the contributions from both 3N forces 
and many-body processes consistently. Our work is based 
on an improved treatment of 3N forces [23.-25] , where the 
normal-ordered one- and two-body parts of 3N forces are 
included to third order in MBPT (in the first calculations 
with 3N forces [22|, their contributions were included 
only to first order) . At the level of particle-particle and 
hole-hole ladder contributions, which corresponds ap- 
proximately to the level of EDF calculations, we find that 
3N forces lead to a reduction in pairing strength compara- 
ble to that seen in the EDF results J^] , although with an 
incorrect odd-even staggering of A„ ' . This inverted stag- 
gering can be related to different mean fields in each ap- 
proach, which has also been observed in Gorkov-Green's 
function calculations [33|- When particle-hole contribu- 
tions are included to third order, reasonable agreement is 
obtained with experimental three-point mass differences. 
Finally, we study the 2+ excitation energies based on the 
same calculations of the even calcium isotopes, and ex- 
tend the discussion of shell evolution to signatures of shell 
closures in the pairing gaps. Our improved treatment of 
3N forces leaves unc hang ed the prediction of A^ = 28 as 
a magic number [23, UM j but reduces the 2+ excitation 
energy to 1.7 — 2.2 MeV in ^'^Ca, which has been recently 
investigated at RIKEN [38,]- 



II. MICROSCOPIC VALENCE-SPACE 
INTERACTIONS AND 3N FORCES 

We calculate the interactions among valence neu- 
trons as in Refs. [23l - [25| . where the theoretical inputs 
to the valence-shell Hamiltonian, single-particle energies 
(SPEs) and residual two-body interactions, are calcu- 
lated microscopically from NN and 3N forces, without 
adjustments. Our results are based on chiral effective 
field theory (EFT) interactioiis, which provide a system- 
atic basis for nuclear forces [33] ■ At the NN level, we 
perform a renormalization-group evolution (4Cl |4l| of the 
500 MeV N^LO potential ^^ to a low-momentum cut- 
off A = 2.0 fm^^ (with smooth re gula tor |43| ) to im- 
prove the many-body convergence |4l| . The resulting 
low-momentum interaction T^ow k is used to calculate the 
two-body interactions among valence neutrons to third 
order in MBPT in a harmonic-oscillator basis of 13 ma- 
jor shells with huj = 11.48 MeV, following the formalism 
of Refs. [13, El. 

At the 3N level, we take into account the leading N^LO 
3N forces [4^ [43] . This includes the dominant long-range 
two-pion-exchange parts (due to A and other excita- 
tions), plus shorter-range one-pion-exchange and contact 
3N interactions. The short-range 3N couplings are fit to 
the binding energy of '^H and the radius of '^He for Viow fe 



TABLE I. Empirical and calculated (MBPT) SPEs in MeV. 



orbital 


empirical 




MBPT 




GXPFl \54} 


KB3G [53] 


P/39/2 


/7/2 


-8.62 




-8.60 


-8.05 


PS/2 


-5.68 




-6.60 


-5.86 


Pl/2 


-4.14 




-4.60 


-3.22 


h/2 


-1.38 




-2.10 


-1.33 


59/2 


(-1.00) 




- 


-1.23 



with A = A3N = 2.0 fm"^ [48|. The fit at lower cutoffs 
approximately includes the effects of the renormalization- 
group evolution in the 3N sector, up to higher-order 3N 
forces. For the two-body interactions among valence neu- 
trons, we include the normal-ordered two-body part of 
3N forces consistently to third order in MBPT. This takes 
into account the interactions of two valence neutrons with 
a nucleon in the core and gives rise to repulsive interac- 
tions between valence neutrons [2l|, [22] ■ These contribu- 
tions are expected to dominate over residual three-body 
interactions based on phase space arguments for normal 
Fermi systems [49| . as shown recently in the context of 
the shell model [50]. The normal-ordered two-body ap- 
proximation has also been benchmarked in ab-initio cal- 
culations ^5l|.[52|. 

The SPEs in ''^Ca are obtained by solving the Dyson 
equation, including one-body contributions to third or- 
der in MBPT in the same space as the two-body interac- 
tions. We include the normal-ordered one-body part of 
3N forces, which corresponds to the interactions of one 
valence neutron with two nucleons in the core. These 3N 
contributions range from 1 — 6 MeV, depending on the or- 
bital, approximately an order of magnitude larger than 
two-body effects, which is consistent with the normal- 
ordering hierarchy [4^, [51J, [S^ • The convergence in the 
one-body sector is slower than the two-body case, but 
the results show convergence in 13 major shells for the 
NN part of nuclear forces, and in 5 major shells for 
the 3N part. At each self-consistency iteration when 
solving the Dyson equation, the unperturbed harmonic- 
oscillator spectrum is updated such that the energy of 
the degenerate valence-space orbitals is set to the cen- 
troid, Ej(2j + l)ej/Ej(2j + 1), of the calculated SPEs 
of the previous iteration. The calculated SPEs are given 
in Table [T| They are comparable to those of the phe- 
nomenological KB3G 53] and GXPFl [5J] interactions 
for the p/-shell orbitals [22]. 

We have previously shown that, in valence-shell cal- 
culations based on nuclear forces, the inclusion of or- 
bitals on top of the standard one major harmonic- 
oscillator shell becomes important for neutron-rich iso- 
topes [22, l2j]. This also is in line with our calcula- 
tional framework, which treats the valence orbitals ex- 
actly by diagonalization, but includes the contributions 
from outside the valence space perturbatively. Therefore, 
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FIG. 1. (Color online) Three-point mass differences A„ from "'"Ca to ^''Ca calculated in the p/-shell to successive orders 
of particle-particle and hole-hole ladders, with and without the leading chiral 3N forces, and in comparison with experiment. 
Panel (a) shows the first-order pairing gaps, while panels (b) and (c) include the effects of adding second- and third-order ladder 
contributions, respectively. All calculations use empirical SPEs in ^^Ca (from GXPFl [5J]). Experimental data are from the 
2003 Atomic Mass Evaluation (AME2003) [gJ and including recent precision mass measurements for ^^■'^^Ca from TITAN [23 ]. 



extended valence spaces treat nonperturbatively as many 
orbitals as possible in the diagonalization. To probe the 
effects on pairing of the correlations involving extended- 
space orbitals, we have performed calculations in both 
the 0/7/2, lp3/2, 0/5/2, lPi/2 valence space (p/-shell) as 
well as the expanded space including the 039/2 orbital 
{pfgg/2 valence space). Both cases start from a ""^Ca 
core. We take two approaches with respect to SPEs: 
in all p/-shell calculations we use the empirical SPEs 
of Ref. [54| . while for the pfgg/2 valence space, we ei- 
ther use these empirical values (with the .gg/2 SPE set 
to — l.OMeV), or the SPEs calculated consistently in 
MBPT. Regarding the pfgg/2 calculations, in TableUwe 
observe that the gg/2 orbital is obtained close to the pf- 
shell, which can lead to substantial contributions from 
this orbital. 

Since the pfgg/2 valence space comprises orbitals be- 
yond one major harmonic-oscillator shell, the center-of- 
mass (cm) motion of the valence nucleons is not guaran- 
teed to factorize in general. As in Refs. [53, Hg], we ex- 
plore possible cm admixture in our calculations by adding 
a cm Hamiltonian, /3Hcm, with /3 = 0.5, to the valence- 
shell Hamiltonian. This results in only modest changes 
of < 300 keV in the pairing gaps. These changes can be 
understood because the non-zero cm two-body matrix el- 
ements are also relevant matrix elements of the MBPT 
calculation, making a separation between these two ef- 
fects difficult. Similarly, we find nonzero (Hem) values, 
which point to possible cm admixture and/or nonneglible 
occupancies of the gg/2 orbital. However, as discussed 
in Ref. [53| and shown in the context of coupled-cluster 
calculations pa\, /3-depcndence and (Hem) ^ do not 
necessarily imply cm contamination. 

Since it is important to understand this issue, work 
is in progress in several directions in both the pfgg/2 
and sdfY/2P3/2 [23| valence spaces. First, we perform 



a nonperturbative Okubo-Lee-Suzuki-Okamoto transfor- 
mation [59l - [63 | to project the extended valence-space 
interaction into the standard major harmonic-oscillator 
shell, which is free of cm spurious states. This will 
nonperturbatively include the contributions from the ex- 
tended orbitals. We will also adapt the singular- value 
decomposition used in Ref. ^58j to test cm factorization 
in the extended-space calculations. Finally, we will ap- 
ply the in-medium similarity renormalization group (IM- 
SRG) 63, 64] to develop extended valence spaces, where 
the cross-shell matrix elements can be suppressed under 
the evolution to the valence-shell Hamiltonian. This will 
ensure (i?cm) ~ and provide a nonperturbative bench- 
mark for MBPT calculations. 

With the developed valence-shell interactions based on 
chiral NN and 3N forces, we perform exact diagonaliza- 
tions in the valence space to obtain the ground states 
and, for the even calcium isotopes, the first excited 2+ 
states. The shell-model codes ANTOINE [13, HI] and 
NATHAN 'm] have been used throughout this work. 



III. RESULTS 

A. Effects of 3N forces and many-body processes 

The EDF pairing calculations of Ref. [la] use the 
same chiral NN and 3N forces of Sect. [Hi where 3N 
forces are included as density-dependent two-body in- 
teractions from Ref. 66], and restricts the contributions 
to pairing to the dominant ^Sq channel. Pairing in the 
EDF approach is realized by solving the Hartree-Fock- 
Bogoliubov equations, which in the pairing channel ap- 
proximately is equivalent to summing ladder diagrams. 
Although the comparison to diagonalization in the va- 
lence space is not one-to-one, it is still instructive to first 
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FIG. 2. (Color online) Three-point mass differences A„ from ''"Ca to ^''Ca calculated to third-order ladders in MBPT with 
empirical SPEs, panel (a), compared with the EDF results of Ref. [isl ]. panel (b). Results with and without the leading chiral 
3N forces are shown following the legend of Fig. [T] and in comparison with experiment [24l.l67l]. 
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FIG. 3. (Color online) Three-point mass differences A„ in the calcium isotopes calculated to third order in MBPT with and 
without the leading chiral 3N forces, and in comparison with experiment [24l.l67l]. The legend is as in Fig. [l] Panel (a) shows 
the results of the third-order ladder contributions. Panels (b) and (c) include all MBPT diagrams to third order in the p/-shell 
and the extended pfgg/2 valence space, respectively. The results in the p/-shell are with empirical SPEs. For the pfgg/2 space, 
we show pairing gaps for both the MBPT and empirical SPEs. 



compare EDF results to valence-shell calculations includ- 
ing particle-particle and hole-hole ladder diagrams only. 

Figure [T] shows the three-point mass differences A„ 
from ''"Ca to ^"^Ca calculated in the p/-shell to successive 
orders of particle-particle and hole-hole ladders. At both 
the NN and NN-I-3N level, increasing the order to which 
ladders are included also increases the pairing gaps, sys- 
tematically improving the agreement with experimen- 
tal data, taken from the 2003 Atomic Mass Evaluation 
(AME2003) [62| and including recent precision mass mea- 
surements for ^I'^^Ca from TITAN [24]. Moreover, the 
results converge rapidly: from first to second order, there 
is a significant increase in A„ ' of ^ 0.4 MeV; from second 
to third order, the change is only ~ 0.1 MeV. 

Pairing gaps calculated at this level are clearly deficient 
with respect to experiment. In addition to being below 
the experimental pairing strength, the odd-even stagger- 



ing of A„ is inverted compared to experiment (A^ is 
stronger for odd masses than for even ones). This in- 
verted staggering is a sign that the mean-field part at 
this level is too attractive, resulting in a lack of satu- 
ration and an incorrect symm etry energy, similar to the 
calculations discussed in Ref. [33j- The correct odd-even 
staggering of A„ ' seen in the EDF results of Ref. [l^ 
(see the right panel of Fig. [2]) is already built into the 
Skyrme functional used in Ref. |15| . 

Taking into account 3N force contributions at the lad- 
ders level, we find in Fig. [1] that the repulsive effect of 

(3) 

chiral 3N forces leads to a systematic suppression of A„ . 
Ranging from 0.2 — 0.5 MeV, this is similar to the de- 
crease in pairing strength observed in the EDF study of 
Ref. [l^ , as can be seen in Fig. [51 Note that the incor- 

(3) 

rect odd-even staggering of A„ persists in the NN-I-3N 
case, but the repulsive 3N forces make it less pronounced. 
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FIG. 4. (Color online) Three-point mass differences A„ in the calcium isotopes calculated to third order in MBPT with and 
without the leading chiral 3N forces, and in comparison with experiment [24l.l67l]. The full calculations (solid lines) are compared 
to the case where all valence neutrons are coupled to J = pairs (seniority s = 0, given by the dashed lines) . Panel (a) shows 
results based on NN forces only in the p/-shell, while panel (b) includes NN-I-3N interactions. Similarly, panels (c) and (d) 
show results in the extended pfgg/2 valence space including NN and NN-I-3N interactions, respectively. All calculations use 
empirical SPEs in ''^Ca (from GXPFl \m 



Moreover, as shown in Fig. [31 the impact of 3N forces on 
pairing gaps varies more along the isotopic chain, when 
all third-order many-body processes are included. 

When particle-hole contributions are included in a full 
third-order calculation, we find in Fig. [3] a clear improve- 
ment compared to including only ladder diagrams. In the 
p/-shell, the three-point mass differences are increased, 
leading to reasonable agreement with experimental data. 
This clearly demonstrates the importance of particle- hole 
many-body processes, such as core-polarization, on pair- 
ing in nuclei. Our results show that they can provide the 
missing pairing strength required to reproduce experi- 
ment on top of the direct NN-I-3N interactions. Analo- 
gously, the systematic differences between theoretical and 
experimental pairing gaps found in the EDF approach of 
Ref. [19] may be attributed to these effects. 

From Fig. [3] (b) and (c) we see that 3N forces have 
a limited effect in the full calculation from '*^~''^Ca 
and ^^~^^Ca. In the first region, experimental three- 
point mass differences are already well described with 
NN forces only. More noticeable differences occur from 
48-54Qg^^ where the addition of 3N forces significantly 

modifies the structure of A„ . In particular the rather 
flat trend in pairing strength from ^^"'^^Ca given by NN 
forces evolves to a structure with peaks at the even iso- 
topes, reproducing qualitativel y t he experimental odd- 
even staggering in this region [2J]. Extending the cal- 
culation from the p/-shell to the pfgg/2 valence space 
provides further improvement when compared to experi- 
ment, in particular with respect to the pronounced peaks 
at ''^Ca and ^^Ca and the corresponding shell closures. 
On the other hand, this overestimates A^ at ^^~^^Ca, 
but most likely this is due to the influence of the some- 
what overpredicted shell closure at '*^Ca rather than a 
deficiency in the pairing channel. We will discuss shell 
structure further in Sect. IIII CI 



B. Seniority truncations 

We can further study the physics of the pairing gaps 
by considering a restricted calculation, in which the full 
third-order wavefunctions are constrained to have all va- 
lence neutrons forming J — Q pairs (seniority s — Q). 
This is studied in Fig. H both for NN and NN-I-3N inter- 
actions in the p/-shell and the pfgg/2 valence space. In 
most cases, at the s = level there is already a reasonable 
description of the three-point mass differences, indicating 
that pairing gaps are, as expected, largely described by 
J = interactions. This is also reflected when using phe- 
nomenological interactions fSJ, [SJI , where the difference 
of the s = to the full calculation is < 100 keV. However, 
for mid-shell isotopes around ^°Ca there is a clear differ- 
ence for the interactions studied here, especially when in- 
cluding 3N forces in the pfgg/2 valence space. This gives 
us an indication how the gg/2 orbit is affecting physics 
in this region, as we find that the main difference in go- 
ing from s = to s = Smax is in the occupancy of this 
orbital. At the level s = 4, where all but 4 particles are 
constrained to form J = pairs, the results are nearly 
indistinguishable from those of the full calculation. 



C. Pairing gaps and shell closures 

A well-known signature of a closed-shell nucleus is a 
relative peak in the three-point mass difference. There- 
fore, this work extends the discussion of the evolu- 
tion of shell structure in the neutron-rich calcium iso- 
topes [22, [2^ • We study the 2+ excitation energies of the 
even calcium isotopes based on the same calculations as 
for the three-point mass differences presented here and 
in Ref. [2J], where 3N force contributions are included 
consistently to third order in MBPT. The resulting 2+ 



>> 3 

u 

w 



rN 



1 7 

- 



-III 


1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 _ 

A ~- 
1 1 1 1 1 1 1 1 1 1 1 


-•— ■ - 
1 • 

, , 1 


■ • 

• 

- NN 

- NN+3N (emp) 

- NN+3N (MBPT) 

1 1 1 1 1 1 1 1 1 1 1 1 1 



44 48 52 56 60 64 68 
Mass Number A 



FIG. 5. (Color online) First 2+ excitation energies in the 
even calcium isotopes with and without 3N forces compared 
with experiment. The energies are calculated to '^^Ca in the 
extended pfgg/2 valence space, using both empirical (dot- 
dashed and dashed lines) and microscopic MBPT SPEs (solid 
line). Experimental energies are taken from Ref. [6^ ]. 



excitation energies are presented in Fig. [5l 

From Fig.[3](b) we observe that the pairing gaps calcu- 
lated in the p/-shell largely follow the results of Ref. [23] , 
where no apparent closure is seen in ''^Ca with NN forces 
and limited enhancement when 3N forces are included. 
Furthermore, with the addition of 3N forces, a very mod- 
est peak emerges at ^^Ca, and an initially strong peak at 
^^Ca becomes strongly suppressed. It is very interesting 
that at ^^Ca, A„ obtained from a combination of masses 
from the recent precision measurements of ^^'^^Ca from 
TITAN nil and the AME2003 extrapolation ^ displays 
a signature of a shell closure at ^^Ca, consistent with the 
high experimental 2+ energy. However, direct mass mea- 
surements beyond ^^Ca are critical. 

In the pfgg/2 extended space calculations of Fig. [3] (c) , 
we find a clear enhancement at "^^Ca with 3N forces, sug- 
gesting that correlations involving the gg/2 orbit are im- 
portant in predicting this shell closure. In Fig. [5l we 
similarly find a high 2+ excitation energy, which is some- 
what overpredicted. The best agreement with experi- 
ment in Fig.[5]is found for the NN-F3N case with MBPT 
SPEs, where all inputs to the valence-shell Hamiltonian 
are calculated based on chiral NN and 3N forces without 
adjustments. In ^'^Ca, Figs. [3] and [S] show that 3N forces 
lead to a quenching of the pairing gap strength and a low- 
ering of the 2+ energy (to 1.7 — 2.2 Me V), regardless of 
the valence space or which SPEs are used. This is in con- 
trast to the calculations of Ref. [22 , based on first-order 
3N force contributions. This shows that the structure of 
A^ = 34 is sensitive to theoretical details of the calcula- 
tion, making this a very good test case, as discussed in 
Ref. [28| . Further investigations are necessary to explore 
the theoretical uncertainties from the truncation of 3N 
forces and to the valence space. 

As we approach ^"Ca, we see no further signatures of 



shell closures in the pairing gaps, except for the interest- 
ing case of ^°Ca, which corresponds to a potential iV = 40 
magic number. In this case, regardless of whether 3N 
forces are taken into account, there is a significant peak 
in A„ . This is supported by the results of Fig. [S] where 
the calculated 2+ energies are approximately twice the 
2+ energies in ^®Ca and ^^Ca. This result should, how- 
ever, be taken with care, because of the theoretical un- 
certainties from the truncation of the valence space and 
from neglected continuum effects, which both are more 
severe for the most neutron-rich calcium isotopes. Re- 
cent coupled-cluster calculations indicate that due to the 
continuum, the higher-lying 1^5/2 and 2si /2 orbitals, not 
considered in this calculation, are lowered significantly at 
^"Ca, both appearing below the 59/2 orbital [28j |. 

Finally, it is worth emphasizing that in Fig. [2l while 
the ^°Ca and to a lesser extent the ^^Ca shell closures are 
evident in the EDF results, there is little indication of 
shell structure throughout the rest of the isotopic chain. 
Indeed, there is also no clear peak in A„ ' for either '"^^Ca 
or ^''Ca in the ladder results of Fig. [51 While the addition 
of 3N forces does strengthen the peak signature in ^^Ca, 
the presence of the peak with NN forces is mainly due 
to the very large spacing of the 77/2 and P3/2 orbitals 
(~ 4MeV) in the EDF SPEs [69|. 



IV. SUMMARY 



We have studied the role of 3N forces and many-body 
processes in a microscopic description of nuclear pairing 
in the calcium isotopes. At the level of particle-particle 
and hole-hole ladder contributions, the repulsion from 3N 
forces decreases the pairing gaps, in line with the EDF 
calculations of Ref. [15|] . Taking into account many-body 
processes to third order in MBPT, we have found rea- 
sonable agreement with experimental three-point mass 
differences. This shows the important role of attractive 
particle-hole contributions to pairing in nuclei [32I435J in 
our first calculation based on chiral NN and 3N forces. 
Finally, we have studied the 2+ excitation energies of 
the even calcium isotopes in the same framework. This 
improves the first calculations with 3N forces [24I, by 
including 3N force contributions consistently to third or- 
der in MBPT as in Ref. [2^. This improved treatment 
leaves unchanged the prediction of iV = 28 as a magic 
number [22, 123, but reduces the 2+ excitation energy to 
1.7 — 2.2 MeV in ^^Ca, which has been recently investi- 
gated at RIKEN [SSJ. 

Although this study focused on the calcium isotopes, 
extensions to heavier isotopic chains such as nickel and 
tin are possible. They will serve as a test of the presented 
valence-shell calculations in heavy nuclei, while providing 
a further link between the shell model and EDF methods 
based on nuclear forces. 
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